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ABSTRACT: Cosmological solutions with a homogeneous Yang-Mills field which oscil- 
lates and passes between topologically distinct vacua are discussed. These solutions are 
used to model the collapsing Bartnik-McKinnon gravitational sphaleron and the associated 
anomalous production of fermions. The Dirac equation is analyzed in these backgrounds. 
It is shown explicitly that a fermion energy level crosses from the negative to positive 
energy spectrum as the gauge field evolves between the topologically distinct vacua. The 
cosmological solutions are also generalized to include an axion field. 



1. Introduction 

A few years ago the Einstein- Yang-Mills field equations were shown to admit spher- 
ically symmetric particlelike solutions in which the Yang-Mills configuration is bound by 
gravity in unstable equilibrium These solutions correspond to gravitational analogs 
of the standard model sphaleron solution 0. In other words, they are saddlepoint so- 
lutions that lie on an energy barrier in field configuration space separating vacua with 
different topological charge. Just as the ordinary sphaleron mediates transitions between 
these vacua and leads to anomalous fermion production so should its gravitational analog. 
Consider a gauge field which vanishes as t — oo and is given by a pure gauge with 
unit Chern-Simons number as t ^ oo. According to the anomaly equation, the change 
in fermion number is given by the change in Chern-Simons number: An = J F A Fd^x. 
Ordinarily, even in a time-dependent background particles and anti-particles are created 
in pairs. However, as discussed in [Q, this is not necessarily the case for fermions. For 
slowly varying fields, the anomalous production of fermions has a simple description in 
terms of the spectral fiow of the Dirac operator. Consider the spectrum in the background 
of the sequence of instantaneous static field configurations with the initial and final gauge 
field configurations being pure gauge but with their Chern-Simons numbers differing by n. 
The initial and final spectra should coincide with the spectrum in the absence of a gauge 
field. However, the anomalous production is realized by n energy levels crossing from neg- 
ative to positive values of the energy. Thus, assuming one has begun in the vacuum with 
the negative energy Dirac sea filled, afterwards, there would be n positive energy levels 
occupied corresponding to the creation of n particles. Since the sphaleron configuration 
is the midpoint in the transition, one would expect that there would be a zero energy 
normalizable fermion mode in the sphaleron background. Indeed, this was shown to be 
the case for the Bartnik-McKinnon solution 0]. 

For the sphaleron of the standard model, one may view the time dependent process 
as being initiated by a high-energy collision of particles and ending with the fields going 
off to infinity. For the gravitational sphaleron, in addition to this scenario, there is the 
possibility of the configuration undergoing gravitational collapse and forming a horizon and 
singularity. Thus, the anomalously produced fermions could either enter the horizon or go 
off to infinity. [One should note that the third possiblity of the fermions remaining bound 
outside the black hole is excluded by the no-hair theorem.] To calculate the fraction 
of fermions that go off to infinity, one should study the evolution of the fermion zero- 
mode. Thus, one should solve the time-dependent Dirac equation in the background of 
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the collapsing Bartnik-McKinnon solution. This is of course a difficult problem for which 
little progress can be made analytically. 

In this paper, we study cosmological solutions to the Einstein- Yang-Mills equations 
with closed spatial slices with the motivation that they might serve as an approximation 
for the interior of the collapsing sphaleron. Since SU{2) and the spatial sections are both 
three-spheres, there is a natural ansatz for a homogeneous SU{2) gauge field involving a 
single function of time, /. The energy momentum tensor of the gauge field corresponds 
to pure radiation, and the spacetime to the Tolman universe. In contrast with abelian 
electromagnetism where the only way to obtain a pure radiation stress tensor is by some 
averaging procedure, with a non-abelian gauge field, one can obtain a pure radiation stress 
tensor directly from a classical solution. As we will discuss, these cosmological solutions 
share features of the sphaleron. The time evolution of / corresponds to motion in a 
double well with minima at the topologically distinct vacua. There is an unstable gauge 
field configuration at the top of the energy barrier which corresponds to the reduction 
of the four-dimensional meron solution to three-dimensions. There is a fermion zero 
mode in the background of this configuration with the same internal structure as the zero 
mode in the Bartnik-McKinnon solution. Moreover, we show explicitly that as the field 
evolves between the topologically distinct vacua there is a single fermion level crossing 
corresponding to the anomalous production of one particle. 

In Section 2, we review Robertson- Walker spacetimes with a homogeneous Yang-Mills 
field. In Section 3, these solutions are extended to include an axion, or antisymmetric 
tensor field. In Section 4, we begin our study of fermions in this homogeneous time- 
dependent background. Since the spatial slices are three spheres (S^), spinor harmonics 
on are reviewed. The coupling to the gauge field is then included. We find a fermion 
zero mode in the background of the / = 1/2 saddlepoint solution. It is shown exlicitly that 
as the gauge field evolves between vacua, a negative energy level arises from the Dirac sea, 
crosses zero energy at the saddlepoint configuration, and then enters the positive energy 
spectrum. 

2. Yang-Mills Cosmologies 

In this section, we consider cosmological solutions with a homogeneous Yang-Mills 
field as energy source 0. Since we wish to match the interior cosmological solution onto 
an exterior solution initially at rest, we take closed spatial sections because only for k = 1 



2 



is there a point at which the solution is momentarily static. Assume a Robertson- Walker 
ansatz 



(2.1) 

= a''{r]){-dTf + dnl) 

with spatial sections where (iOl is the round metric on 5""^ and where rj conformal time 
with ^ = . The metric (|2.1| ) is conformal to the Einstein static universe. In addition to 
the spacetime metric, consider an SU{2) Yang-Mills field where a is a Lie algebra index 
with associated field strength F = dA — ze/2[A, A] where A = A^dx^ is the Lie-algebra 
valued one- form and e the gauge coupling constant. Under a gauge transformation, U{x), 
the fields transform as A — U~^AU i/e U~^dU and F — > U~^FU. As we now discuss, 
since the group manifold of SU{2) is and coincides with the spatial sections, there is a 
natural ansatz for the gauge field which shares the symmetries of dQ^, the spatial metric. 

Before discussing the gauge field ansatz, we review invariant forms on S^. The natural 
imbedding of into the group manifold of SU{2) is given by 

9=[f/X' f^-%)=^"+'^'^'^ detg = {x'r+{^'rHxrH^r = i (2.2) 

where r* are the Pauli spin matrices. 5*0(4) transformations are induced by its two fold 
cover SU (2) l x SU (2)fj corresponding to multiplication of g on the left and right by SU (2). 
Left invariant one-forms are given by (or more precisely, given by the pull back to of) 

ef = -iTrT^g-^dg = 2{x^dx' -x'dx^ + e^JkX^dx^), i = l,2,3 (2.3) 

and obey def = ^eijkej A ef. Under right multiplication the forms ef (|2.3| ) transform in 
the adjoint. The right invariant one-forms are ef^ = iTr r'^dgg~^ and obey de^ = ^eijkef /\ 
ef-. ef and ef can be written as ef = eijk{^jk + *^jk) and ef = eijk(M.jk — *Mjfc), the 
self-dual and anti-self-dual parts of M^j = x^dx^ — x^dx'^ where *Mij = ^tijki^ki- The 
left and right invariant vector fields are dual to ( p.3|) and given by 



(2.4) 



^f = V2(.^|,-x^|, + e,..^|,) 



dx^ dx'^ dx'' 



Ef and generate the transformation {z^,z'^) (e^^/^z"*^, e*^/^z^) and {z^,z^) —>■ 
^g-i6i/2^i^ respectively where = x'^ + ix^ and z'^ = x^ + ix"^. Li = —iEf and 

Ri = —iEf^ obey the SU{2) commutation relations [Li,Lj] = ieijkLk, [Ri,Rj] = ieijkRk, 
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and [Li,Rj\ = 0. The round metric is the puU back to of the flat four dimensional 
metric 

ds^ = {dx^f + {dx^f + {dx'^f + {dx^f (2.5) 



and in terms of the matrix g (|2.2|) takes the form 



ds'^ = -ilV {g-^dgf = ^ef ® ef = ^ef ® ef . (2.6) 

( ^^) is the bi-invariant metric with isometry group SU{2)l x SU{2)ii corresponding to 
left and right translations. 

There is a natural ansatz for the gauge field on 5''^ which takes the form 

A=-fiv)g-'dg = --fiv)e^^. (2.7) 

Since (|2.7D is invariant under left translation and changes by a global gauge transformation 
under right translation, the gauge field is homogeneous sharing the symmetry group of the 
5""^ spatial slices. / and rj are dimensionless while has units of [M]~^[L]~^. The field 
strength associated with ( |277| ) is 



F^^, F^ = -^dr7Ae^^ + l/(/-l)eV,e^^Ae^^ (2.^ 



Prom this expression or directly from ( |2.7| ), one observes that / = and 1 correspond to 
vacuum, or pure gauge field, configurations. For a general (non-vacuum) configuration, 
the Chern- Simons number is given by the integral Ncs = /gs ^^3 of the three- form 

cj^= Tr{AAdA-—AAAAA) (2.9) 
8tt^ 3 



satisfying doos = -^TrF A F. Substituting (|2.7|) in (|2.9| ) and using the fact that (^ 



implies that | / ei A 62 A 63 = 2tv^, the area of the unit 3-sphere, one eventually finds 
Ncs = 3/^(1 — 2/3 f). The / = 1 configuration has unit topological, or Chern-Simons, 
number corresponding to the fact that the map g has unit winding number. 

The gauge field should satisfy the Yang-Mills field equation *D * F = D^F^^ = 
where D = d—ie[A, ] is the gauge covariant derivative and * is the Hodge dual. Substituting 
in the ansatz (|2.7|) , one obtains an equation of motion for / 



drf df 



— y, V = 2{r-fr. (2.10) 



The field equation is independent of tlie scale factor a{r]) since the Yang-Mills field equation 
is conformally invariant. ( p.lO| ) describes a particle moving in a double well potential V 
with minima at the pure gauge configurations / = and 1. These are the only pure 
gauge configurations within this ansatz. There is a conserved (dimensionless) energy type 
quantity E = ^/^ + V{f) for the equation of motion ( |2.1C1| ). From the form of the field 
strength above, one can identify the electric and magnetic field strength with the kinetic 
and potential energy respectively. In addition to the pure gauge configurations, the / = 1/2 
configuration at the top of the energy barrier with E = Eq = 1/8 is also a static solution 
which is unstable and has non-zero (purely magnetic) field strength 0. This solution in 
fact corresponds to the reduction to three dimensions of the four dimensional Euclidean 
Yang-Mills meron solution |^]. Using the conserved energy, E, one can integrate to obtain 
/(?]) in terms of elliptic functions. For E = Eq, there are two solutions starting and 
finishing at the top of the energy barrier that take the simple form 

The coupling of the Yang-Mills field to gravity is governed by the energy-momentum 
tensor 

^ 1 

= E(^Ma^r - -^g.^K^F^''^)- (2.12) 

1=1 

Thus, each term in ( |2.12| ) behaves locally like the stress tensor of a Faraday fiux tube with 
longitudinal tension and perpendicular pressure. For our ansatz, the three fiux tubes are 
mutually orthogonal and have equal strengths. Thus, the energy-momentum tensor is of the 
form of pure radiation P = p/3 where the energy density is given by p = po/a'^, po = 12^. 
The solution to Einstein's equation for pure radiation is a Tolman universe with scale factor 
in conformal time given by a{r]) = {327rEy/% sinry. Here, Iq = G^l'^e ^ defines the only 
length scale in the system. By contrast, mo ~ G~^l'^e~^ sets the scale of the mass of the 
Bartnik-McKinnon solutions. 

The scale factor aij]) has the familiar a ~ t^l'^ form for small time for a radiation 
dominated universe. These solutions describe an expanding universe of radiation in whcih 
the Yang-Mills field constituting the radiation oscillates in the potential well given in ( |2.10|) . 
We note that even though / = 1/2 is a static solution to the Yang-Mills equation with the 
field strength, F^,^, independent of time, the scale factor is time dependent. The lifetime 



of the universe is A?] = tt, or At ~ \fElo- The number of osciUations in the potential weU 
in the hfetime of the universe can be obtained by integrating /. For large one has 

A77 = / , oc E-^'\ (2.13) 

/ ^2{E-V{f)) 

the integral being taken between the turning points. Thus, the number of oscillations in 
the lifetime of the universe grows as ~ E^^^ for large E. As E tends to Eq from above, 
the period of an oscillation diverges as A77 [E — Eq)~^/'^. 

As stated earlier, one motivation for studying these solutions is that they might serve 
as a model for the interior of the collapsing Bartnik-McKinnon sphaleron just as a simple 
model of the gravitational collapse of a star is provided by a portion of a /c = 1 pressure free 
Priedmann universe. There are, however, certain difficulties with modelling the interior by 
a homogeneous Yang-Mills field. First, unlike the case of a pressure free perfect fiuid, even 
if the stress tensor were initially homogeneous and isotropic, there is no reason to suppose 
that it will remain so. Second, the Bartnik-McKinnon solution more closely approximates 
a shell of Yang-Mills gauge field rather than a ball, and numerical studies suggest that 
the shell becomes thinner as the configuration collapses. Third, since the Yang-Mills field 
possesses non-zero pressure, there will be a pressure jump, and therefore, unlike pressureless 
dust one cannot simply attach the homogeneous interior directly onto the exterior vacuum. 
This last problem might be resolved provided one could smooth out the boundary between 
the homogeneous solution and the exterior vacuum. 



3. Axion Field 

In this section, we consider the addition of an axion field, H^^p. The natural ansatz 
for this field which shares the symmetries of is given by 

H = h{r])e^ (3.1) 

where is the volume form on with metric dQ^. The equations of motion for H are 



d* H = and dH = 0. (3.1) automatically satisfies the first, while the second implies that 



/i is a constant. The energy-momentum tensor for H is 

T^iy = H^xpHj^^ - -Qfj^i.HaiSjH"^"' . (3.2) 

Substituting in (|3.1|), one finds that T^j^ is of the form of a perfect fluid with P = p = h'^/a^. 
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Let us now find the solution for the scale factor, a, with the combined Yang-Mills and 
axion fields as energy source. The Friedmann equation is 



or in conformal time 



a J 1 SttG 

+ 72 = (3-3) 



da ^2 2 SttG 4 



where the energy density, p, is the sum of the Yang-Mills and axion contributions p = 
Pym + Paxion = 12E / a~'^ + h'^a~^ . (|3.4[) can be integrated exactly to yield the scale 



factor in closed form 

1/2 



a(7?) = ((^2 +/3)V2gin2(r7-77o) + a) 



. „ a IQnGE ^ SttG , o (3.5) 

sm 2770 = — a= — j3 = —-h , . 

_^ ^ 3 

The solution still has an expansion followed by a contraction. Near the singularity, however, 
the contribution of H to the energy density dominates with the scale factor vanishing as 
a ~ t^/^ rather than a ~ t^/^ as in the Tolman universe. 



4. Fermions 

We now consider fermions propagating in the Einstein- Yang-Mills backgrounds dis- 
cussed above . The Dirac equation for a massless Dirac fermion coupled to gravity and 
a Yang- Mills field is given by 

«7^(V^ - ieAlT'')^^ = (4.1) 

where is the covariant derivative and are generators for the gauge group. The 
massless Dirac equation is conformally invariant. That is, given a solution , g^,^, A^) in 
d spacetime dimensions, there is another solution 0~2(7^,^, A^). Thus, if \E' is a 

solution in the spacetime ( |2.1|) , ^ = a'^/2\l/ is a solution to the Dirac equation on 5*"^ x R, 



the Einstein static universe (ESU). The Dirac equation (|4.1| ) on ESU takes the form 



(^7°|- + ^^3 +e4"T«)^' = (4.2) 
orj 
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where Ip^ is the Dirac operator on S^. Consider the chiral representation for the gamma 
matrices 

^0 1\ , [ Q 5 f-l 



where cr* are Pauh spin matrices. 

In the absence of the gauge field, solving reduces to finding spinor harmonics on 
5""^. Before considering spinors, we review scalar harmonics on S^. The scalar spherical 
harmonics satisfying V^fi = Xifi are given by 

fi = ai-^...i^n'^^ ■ ■ ■ = x^/r, ct^a^ = (traceless) 

(4.4) 

Xi = -1(1 + 2) 

where V3 is the Laplacian on S^. This is easily checked using the form of the four- 
dimensional flat space Laplacian in spherical coordinates 

^2 3 a 1 2 

'^' = d^+rd~r'^^^' (4-5) 
and using V|(r^/;) = from ( [4.4|) . [In n-dimensions, the harmonics are also given by 
(g), but with eigenvalues A; = -/(/ + n- 1).] Since SU{2)l x SU{2)r is the two-fold 
cover of 5'0(4), cti^.-.i, can be written in terms of spinor components as cx-a^Ax-'-AiAi where 
undotted and dotted components correspond to left and right factors of SU (2) respectively. 
The fact that CKi^.-.i; is totally symmetric and traceless implies that (y-A-i^A^---AiAi '^^ totally 
symmetric in dotted and undotted indices separately. Since the | representation of SU (2) 
is given by rank k totally symmetric SU{2) tensors, the harmonics transform in the (^, |) 
representation of SU{2)l x SU{2)r. [Note that the harmonics for a squashed sphere form 
a representation of some subgroup of SU{2)l x SU{2)r.] Prom (|2.4|) , one can show 

E^E^ = E^E^ = Vl (4.6) 

If one diagonalizes L3, and with eigenvalues m and n, then the harmonics can be 
expressed as monomials in 2:1,2:2,^1,^2 of degree /. The eigenvalues m and n of these 
monomials are determined from those of 21,22,21,22- They are m = and 

10 2 ' 2 ' 2 ' 2 ■ 

Having constructed the scalar harmonics, we now consider spinor harmonics on 5""^. 
Like the scalar harmonics, spinor harmonics are most simply expressed in terms of quan- 
tities in the four-dimensional imbedding space. The four-dimensional Euclidean flat space 
Dirac operator in spherical polar coordinates is 

^ =7'^(|- + |-) + -^3 . (4.7) 
or 2 r r 



The spinor harmonics satisfying Ip^ = are given by 



A. = -(/ + ^) 



where ^T^ - i; constant with a a Dirac spinor index and i a vector index. This is easily 
checked using ( F^ ) and that (0) imphes ^ [r^-^i) =0. T^^.. ■ can be decomposed into 



tensors T^AiAi---AiAi with / + 1 undotted and / dotted indices and tensors T^AAi^i-.-Aiyii 



with / undotted and / + 1 dotted indices. (^4.8|) imphes that they are totaUy symmetric 
tensors. Therefore, the spinor harmonics faU into (| + |, |) and (|, | + ^) representations 
of SU{2)l X SU{2)ii with / = 0, 1, ■ ■ -. Since one has {^3 ,7*^} = 0, the harmonics with 
positive eigenvalues can be obtained by applying 7'^ to ( |4.8D . 

The spinor harmonics can also be expressed in terms of an orthonormal frame on 
S^. The Dirac operator is then given by 

tiPs =^7^^3^, Ds^ = d^ + ^u;'^E^\ E^^' = -[^ , 7^]/4 (4.9) 

where E*-' are the Lorentz generators in the spinor representation with the Dirac matrices 
( ^?3|) satisfying {7*57-'} = —26^^ and u^J , the spin connection. As usual Latin and Greek 
indices refer to tangent and curved space and are related by e^. If we let the left-invariant 
vector fields ( p.4| ) = 2E^^ define the orthonormal frame, then, from (2^), the dual one- 



forms = ^e^* satisfy dei = eijk^j A e^. We can now read off the spin connection Uij 



eijk^^ ■ Substituting the spin connection and the Dirac matrices in ([4.9| ), one eventually 
finds 

-i^3 =(_^Q ^°)' Ho=^S-L + ?,/2, S-L~S'V (4.10) 

with 5"* = cT*/2 and L* = —iE\^. Hq commutes with Ri and Ji = Li + Si. The spinor 
harmonics fall into the (| + ^, |), / = 0, 1, ■ ■ ■ and (| — |, |), Z = 1, 2 ■ ■ ■ representation of 
SU{2)j X SU{2)ji and can be labelled by (/, n,j± = | ± l^is)? the eigenvalues of R^, i?3, 
, and J3. They take the form of products of 2-component spinors and scalar harmonics. 
Since S ■ L = {J^ - L'^ - S^)/2 = + 1) - |(f + 1) - |)/2, we find that the spectrum 



2 V 2 

of is 



I . 1 



Ho\l,nJ±j3 >= /if K,n, j±, >, i± = 2 =^ 2 

/if = / + | / = 0,1... (4.11) 
//- = -/-^, / = 1,2,.... 

9 



We thus have recovered the spectrum ( [4.8|) . The spectrum of Hq is symmetric about 
zero. The lowest lying states in the positive energy spectrum are |0,0, |,±| >= |± > 
(S'sli >= ^|± >) with eigenvalue | corresponding to the two components of the constant 
spinor (with respect to the left invariant frame) . The highest energy states in the negative 
energy spectrum are ||,^,0, >= (^i|+ > —Z2\— >)/\/2 and |^,— 1,0, >= (^2!+ > 
—zi\— >)/ \/2 with eigenvalue — | . 

We now consider the Dirac equation ( [4 .21 ) with a non-zero gauge field background. 
Let the fermions transform in the fundamental representation of the gauge group so that 
rpa _ 12 ^ii\y q-O'^ i\yQ Paull splu matrices. From (|2.7| ), we have = — 2//e5* implying 
cT*ef A"T" = —Af/eS ■ T. The full time-dependent Dirac equation (|4.2|) now becomes 



^^ = T^^±, H = Ho + Af{r^)S-T (4.12) 

with ^0 given in (|4.10| ) and where we have decomposed ^' = ^ ^ into i/'-i- , left and 
right chiral components each carrying two-component Lorentz and internal indices. For a 
vanishing gauge field {f = 0), H reduces to Hq. The / = 1 configuration is pure gauge 
with Chern-Simons number unity, and thus, the spectrum of H should be identical to 
Hq. Indeed, we can redefine L' = L + T, and the spectrum in the absence of the gauge 
field is recovered. Since Hq and the interaction Hi = 4fS ■ T do not commute, one 
cannot in general solve ( [4.12| ). However, it is possible to do so within the s-wave sector 
where Li = Ri = 0. In this sector, there are the singlet (or hedgehog) xi and triplet X2 
states satisfying (5* -|- T^xi = {S + T)'^X2 = "2x2- The time-dependent single and 
triplet wave functions are found from ( |4.1UD and ( [4. 121 ) and by performing the conformal 
transformation back to the Robertson- Walker spacetime |^]. One finds 

V'i^ = a-'/Vexp±z f Ei{r,)dr], Ei{f]) = ^ - 3f{r]) 

r 3 ^^-''^ 

'^t^a ^/\2exp±z / ^2(77)^77, ^2(7?) = - + f{ri) 

Let us now consider the adiabatic approximation where the Hamiltonian, H, is treated 
instantaneously. The energies of the single and triplet states are given by Eiirf) and E2{rj). 
Now as / evolves between the vacua / = and / = 1, we see that the singlet (hedgehog) 
state passes from the lowest lying state in the positive energy spectrum, i?i = |, to the 
highest lying state in the negative energy spectrum, Ei = — |. At the midpoint of the 
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transition corresponding to the unstable / = 1/2 solution, the hedgehog xi is a zero 



energy mode [10|. As pointed out earlier, / = 1/2 is the analog of the sphaleron, and like 
the Bartnik-McKinnon sphaleron, it has a fermion zero mode. Since E2 remains positive 
for < / < 1, there is no level crossing in the triplet state. In fact, since the interaction 
25" ■ T is bounded between —3/2 and 1/2, among all sectors the only level crossing occurs 
for the singlet state. This is consistent with the anomaly equation according to which if the 
Chern-Simons number of the gauge field changes by unity, there will be one level crossing 
in each chiral sector. For Dirac fermions, the level crossing cancels between chiral sectors. 

We now consider the effect of the antisymmetric tensor field on the fermions. Consider 
an interaction of the form 

S = ibj ^H^^pri^Y'^Vgd^'x (4.14) 
with b having units [L]^/^[M]~^/^. The Dirac equation now becomes 

(n^ + e^y^'A^^T" + ibH^^p^i^YY) ^ = 0. (4.15) 

One can show generally that a co-closed axion field (H = >Kd^) can be absorbed by a local 
chiral transformation on ^' of the form 

exp (-6z67^(/))^'. (4.16) 

For the axion field, (|3.1| ), one finds that (p = J and therefore, from ( |4.13| ) and ( [4.16| ), 
one finds the time-dependent triplet and single wave functions 

^+ = a-3/\iexp±z [{Ei{r]) - ^)dri 

i^t = a"^/^X2exp±z / (^2(^7) - -^)dV- 

In this paper, we have studied cosmological solutions with a Yang-Mills field as a po- 
tential analytic model of the gravitational collapse of a Bartnik-McKinnon sphaleron. This 
provides a qualitative picture of anomalous fermion production and gives some support 
to the viewpoint proposed in the previous paper Q. However, it clearly cannot give a 
detailed quantitative description. This can, presumably, only be achieved by a numerical 
analysis. 
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